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We propose a class of projected BCS wave functions and derive their parent spin Hamiltonians. 
These wave functions can be formulated as infinite Matrix Product States constructed by chiral cor- 
relators of Majorana fermions. In ID, the spin Hamiltonians can be viewed as SO(n) generalizations 
of Haldane-Shastry models. We numerically compute the spin-spin correlation functions and Renyi 
entropies for n = 5 and 6. Together with the results for n = 3 and 4, we conclude that these states 
are critical and their low-energy effective theory is the SO(n)i Wess-Zumino-Witten model. In 2D, 
we show that the projected BCS states are chiral spin liquids, which support non-Abelian anyons 
for odd n and Abelian anyons for even n. 



PACS numbers: 75.10.Pq, 11.25.Hf, 03.65.Fd 

Introduction. - Efficient description of quantum many- 
body systems is a challenging problem in modern physics, 
as the dimension of the Hilbert space scales exponen- 
tially with the number of particles. For strongly inter- 
acting many-body systems, much of our understanding 
of their properties comes from physically motivated trial 
wave functions and/or exact solutions of specific mod- 
els. A great success of the trial wave function approach 
is the celebrated Laughlin wave function for the frac- 
tional quantum Hall effect at \/m (with m odd) filling 
[ij. Toward exact results, Bethe's solution of the spin- 
1/2 Heisenberg chain fij and integrability of the spin-1/2 
Haldanc-Shastry model [3j] provide invaluable insight for 
critical spin-1/2 chains. 

The justification of trial wave functions is usually a dif- 
ficult task. For example, the relevance of Anderson's res- 
onating valence bond (RVB) state Q for the mechanism 
of high-Tc superconductivity is still a controversial issue. 
A useful technique for justifying trial wave functions is 
to study their parent Hamiltonians for which the trial 
wave functions are exact ground states. For Laughlin 
wave function, the parent Hamiltonian which consists of 
certain Haldane pseudopotentials Q differs from physical 
Coulomb interactions but their difference can be viewed 
as a perturbation [6]. Similar situation arises for the 
spin-1 AKLT state and its parent Hamiltonian [7], which 
contains an extra biquadratic term apart from Heisen- 
berg interactions. Since the spin-1 AKLT model can be 
adiabatically connected to the spin-1 Heisenberg chain 
without closing the gap, it is widely believed that the 
AKLT state qualitatively captures the physics of spin-1 
Heisenberg chain. 

In this work, we propose a class of projected BCS 
states and derive their parent Hamiltonians. These states 
can also be represented as infinite Matrix Product States 
(MPS) Q constructed from chiral correlators of Majo- 
rana fermions. In ID, the spin Hamiltonians are SO(ri) 
generalizations of Haldane-Shastry models. We numer- 
ically calculate the spin-spin correlation functions and 
the Renyi entropies for n = 5 and 6 and compare the nu- 



merical results with field theory predictions from SO(n)i 
criticality. Together with the known results for n — 3 and 
4, we expect that for general n these states are critical 
and belong to the SO(n)i WZW universality class. We 
also show that the projected BCS states with modified 
Cooper pair wave functions provide a good description 
for Ising ordered and disordered phases close to SO(n)i 
criticality. In 2D, the projected BCS states are chiral spin 
liquids with p-\-ip pairing symmetry. We find that these 
topological states support non-Abelian Ising anyons for 
odd n and Abelian anyons for even n, respectively. 

Projected BCS wave function.- Constructing the pro- 
jected BCS wave functions relies on a slave-particle repre- 
sentation of the SO(n) algebra. Let us start from a ID pe- 
riodic chain with even N sites, where the n vectors in each 
site are represented by using singly occupied fermions, 
= c'llO) (a = l,...,n). In terms of fermions, the 
SO(n) generators are written as L"*" — i{c\cb — clca), 
where I < a < b < n. To remove unphysical states in this 
fermionic representation, a single-occupancy constraint is 
required, X]a=i a'^j,a = 1 = 1, . . . , A^, which defines 
a Gutzwiller projector Pq. Then, the projected BCS 
wave function of our interest is defined by 



I*) =PGexp K]. 
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where ^2=i'4a'^]a creates an SO(n) singlet between 
sites i and j. Note that I^P) is a coherent superposition of 
valence-bond singlets of arbitrary range (See Fig. [IJ and 
hence can be viewed as an RVB state If we choose 
Zj = exp(z^j), the amplitude of the Cooper pair wave 
function l/\zi — Zj\ is the inverse of the chord distance 
between the sites. Under this choice, |^) is both real and 
translationally invariant, which is the uniform case that 
we will consider in the following. 

Before discussing the properties of \'^) for general n, 
we establish the relation between ([IJ and some known 
results. For n = 3, after switching to the standard spin-1 



basis = ^(1 - 1) - |1)), \n^) 
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FIG. 1: (color online). Schematic of a valence bond config- 
uration in the projected BCS state ([T]). The valence bonds 
(blue) are SO(n) singles formed by two SO(n) vectors. In the 
uniform case Zj = exp(i22lj)^ the periodic chain is viewed as 
a unit circle and \zi — Zj\ is the chord distance between two 
sites. 

|n^) = |0), we find that the projected BCS state ([T]) 
is equivalent to the spin-1 Haldane-Shastry state, which 
has been considered in Ref. ■ It was shown ^ [13] 

that this state is critical and its low-energy effective the- 
ory is an SU(2)2 (or equivalently S0(3)i) WZW model. 
In a recent work a projected BCS wave function 
similar to ([ij is used as a variational ansatz to describe 
the phases in spin-1 bilinear-biquadratic chain, including 
the Takhtajan-Babujian model |13J which also belongs to 
SU(2)2 WZW universality class [ijj. For n — A, after rep- 
resenting the four S0(4) vectors by two spin-1/2 states, 
we find that the projected BCS state ([T]) can be rewritten 
as a product of two decoupled spin-1/2 Haldane-Shastry 
states. An immediate consequence of this decomposition 
is that the S0(4) state is critical and represents the fixed 
point of the S0(4)i WZW model. 

Infinite MPS and parent Hamiltonian.- From the 
known results for n — S and 4, one may speculate that 
for general n the projected BCS state ^ belongs to the 
SO(n)i WZW universality class. Let us further uncover 
this relationship by formulating ([T]) as an infinite MPS. 
The SO(n)i WZW model has a primary field with con- 
formal weight hy — 1/2 in the vector representation [l5| , 
which can be naturally interpreted as Majorana fermion 
fields x"(z) (a = 1, . . . , n). This Majorana representation 
of the primary field allows us to rewrite the projected 
BCS state ^ as the following infinite MPS: 

n 

|*)= •■,ajv)|n"\...,n""), (2) 

ai ,...,ajv— 1 

where the coefficients are the chiral correlators of N Ma- 
jorana fields [16|] 

vl/(ai, . . . , aAT) = ix""' (^i)x"^ (^2) • • • X"" (zn))- (3) 

Unlike usual MPS with finite matrix dimensions, the 
state is an infinite MPS iS] since its ancillary space on 



which the Majorana fields act is an infinite-dimensional 
Hilbert space. This allows the infinite MPS ^ to de- 
scribe the expected SO(n)i criticality with unbounded 
increase of the entanglement entropy. 

The key benefit of the infinite MPS formulation is that 
a parent Hamiltonian can be derived, such that ([3]) is its 
exact ground state. As shown in Q , the presence of null 
vectors in conformal field theories (CFT) leads to a set 
of operators which annihilate the infinite MPS. Following 
this approach, we have derived such operators for ([3|) 

where my ee {zi + Zj)/{z,-Zj) and Li-Lj = Yja<b^f^f- 
Since Af 1^-) = Vi,a,6 and E^^f I*) = Va, 6, we 
can define a parent Hamiltonian H = E - a<b(A"'')^A"'' -t- 
2(JV-2) Y,a<bmiLff + Ea whose ground state is the 
infinite MPS ^ with energy Eq. Choosing — — |(n — 
1)A^(A^^ — 4), the explicit form of H is given by 

~ 3fa-'l) ^ WijW^kiL^-Lj){L,-Lk). (4) 

Generically, the Hamiltonian ([5]) is a long-ranged SO(n) 
bilinear-biquadratic model with three-spin interactions. 
For n = 4, as we expected, the Hamiltonian only has 
inverse-square Heisenberg exchange interactions, which 
can be decomposed into two spin-1/2 Haldane-Shastry 
Hamiltonians due to SO(4)~SU(2)xSU(2). 

Jastrow versus Pfajfian.- It is well-known that the 
SO(n) Lie algebra has a sharp difference between even 
and odd n [17|. As we shall see, this leads to distinct 
forms of the wave function ([T]) in Cartan basis for even 
and odd n: The former has a pure Jastrow form, while 
the latter includes a Pfaffian factor. To sec this differ- 
ence, let us consider S0(2/) and S0(2^ + 1) (/: integer) 
and choose mutually commuting Cartan generators as 
L^'^,L^^, . . .i2i-i,2i_ Pq^. sO(20, a convenient choice of 
the Cartan basis is defined by |0, . . . , nia — ±1, . . . , 0) — 
dn^") ± i|n2"-i))/v/2 (a = 1,...,0, where m„ is the 
eigenvalue of 7^20-1,20^ p^^. |-j^g vectors \0,...,ma = 
±1, . . . , 0), we label their positions in the spin chain by 
x^^"-* < • • • < xj^''. In this basis, the wave function ([T]) for 
even n — 21 takes the form 

^({m}) = Prn n n(^» - , (5) 

u I (1) (1) (0 (0^ / 

where p,„ = sgn(a;i ' . . . xy^ . . .x\' . . .xy^) (sgn: signa- 
ture of a permutation) if ma^i = Va and pm = 
otherwise. 



3 



For S0(2Z + 1), apart from the vectors |0, . . . , ttiq 
±1, . . . , 0), there exists an additional vector |0, . . . , 0) 



^2l + l\ 



which is annihilated by all Cartan generators. 



Labeling their positions by x^^-* < 
function ([1} for odd n = 21 + 1 is written as 



< x^^\ the wave 
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n n(-' 
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(6) 



where = sgn(xi ' . . .x'j^'^,x\ ' . . .xy^ . . .x\' . . .xy^) \i 
X)i "^a.i — Va and pm — otherwise, and the Pfafhan 
factor Pfo( ^.l^. ) is restricted to the positions for the 
extra vector |0, . . . , 0). 

Numerical results.- The power-law decaying correla- 
tion functions and the universal scaling of entanglement 
entropy 18| are characteristic behaviors of conformal 
critical points in ID. Even though these quantities are 
difficult to be computed analytically for ([T|), the Jastrow 
and Pfaffian forms ([6]) and ([7]) of the wave functions are 
very suitable for determining them numerically via the 
Metropolis Monte Carlo (MC) method Below we 

focus on the projected BCS state ([T]) with n — 5 and 
6 and provide numerical evidence that they belong to 
S0(5)i and S0(6)i WZW models, respectively. 

For critical spin chains in SO(n)i WZW universality 
class, field theory predicts that for n < 8 the spin-spin 
correlation function behaves as ~ {—l)^/x^ 

with 77 = n/4 [i^. For the projected BCS state ^ with 
n = 5 and 6, we have computed the two-point spin corre- 
lator (ij^Lj^^). Fig. [2] shows the numerical results for 
N = 200. The critical exponents that best fit with our 
numerical data are 77 = 1.22 for S0(5) and 77 = 1.42 for 
S0(6) (solid lines in Fig. [5]), which agree very well with 
the field theory predictions (dotted lines). 

The entanglement entropy that is easily accessible via 

MC method is the Renyi entropy S'l '' = — lnTrp| (See 

Ref. [8, 21 23]), where pl is the reduced density matrix 

of the state in a subsystem of length L. For SO(n)i WZW 

(2) 

model with c — n/2we expect S)^ — cln[sin(7ri/Af)]/4+ 
C2 [l8| . where c'2 is a constant. For n = 5 and 6, we plot 
Sf^ as a function of ln[sin(7ri/iV)]/4 for N = 200 in 
Fig. [21 From our MC data, the estimates of the central 
charge are c = 2.31 for S0(5) and c = 2.76 for S0(6) 
(solid lines in Fig. which are close to the predicted 
c = n/2 (dotted lines) but show some deviations. 

The origin of the small deviations of the numerical re- 
sults and the SO(?^)i predictions may be due to the pres- 
ence of marginally irrelevant terms in the SO(n)i WZW 
model for ([T]) and its parent Hamiltonian ([5|), unlike the 



„(0 



SU(n) Haldane-Shastry models [2J, |25[ (including the 
spin-1/2 Haldane-Shastry model for n = 2) being the 
fixed points of the SU(n)i WZW model. For n = 3, the 
presence of marginal term in the spin-1 Haldane-Shastry 
model has been confirmed numerically I, [1^. If this 
is also the case for n > 5, an interesting open question 
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FIG. 2: (color online). 

(logarithmic scale) \n[{ — l)^{Lj'^Lj'^^)] as a function of 
ln[sin(7ra;/iV)] in the projected BCS state gl) for iV = 200 
and (a.) n — 5 and (b) n = 6. The sohd lines (red) are fits 
of the form ln[(-l)^ (Lf Lj^^)] = ri\n[sm{nx/N)] + A, where 
ri and A are fitting parameters. The dotted lines are also fits 
of this formula but with the field theory prediction rj — n/4: 
(Ref. M 
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FIG. 3; (color online). Renyi entropy S^^^ as a function of 
ln[sin(7ri/iV)]/4 in the projected BCS state ^ ioT N ^ 200 
and (a) n = 5 and (b) n — 6. The solid lines (red) are fits 
of the form 5*^^' — cln[sin(7rL/Ai')]/4 -I- c'2, where c and c'2 
are fitting parameters. The dotted lines are also fits of this 
formula but the central charge c is fixed to c = n/2 of the 
SO(n)i WZW model. 



is whether there exist a modified version of ([T]) and its 
parent Hamiltonian that represent the fixed point of the 
SO(n)i WZW model. 

Away from SO(n)i criticality.- After showing that 
the projected BCS state ([T]) captures the physics of the 
SO(n)i WZW model, it is natural to ask whether simi- 
lar projected wave functions are relevant for gapped spin 
chains away from (but close to) SO(n)i criticality. Let 
us restrict ourselves to SO(n) symmetric models for sim- 
plicity. According to the well-known result by Witten 
[26|, the SO(n)i WZW model is equivalent to n mass- 
less Majorana fermions, i.e., n Ising models at criticality. 
For this critical theory, the only relevant perturbation al- 
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lowed by SO(n) symmetry is the mass term of Majorana 
fermions. Thus, the low-energy effective theory has the 
following Hamiltonian density: 



'^T.^n^L: (7) 



i/=i 



where C^(^) are right (left) moving Majorana fermions, v 
and m are their velocity and mass. Here we have assumed 
four-fermion interactions are weak and can be neglected, 
since they are marginal and only renormalize the mass of 
Majorana fermions at low-energy limit 

The SO{n)i criticality corresponds to to = 0. The 
two gapped phases adjacent to the SO(ri)i criticality 
are i) Ising ordered phase (to < 0) and ii) Ising disor- 
dered phase (to > 0). For these two phases, we note 
that they can be well described by modified projected 
BCS states. Actually, these two gapped phases and an 
SO(n)i critical point (Reshetikhin model 28 [) are real- 
ized in the SO(n) bilinear-biquadratic chain 2^, 2^. The 
ideal example that belongs to the Ising ordered phase 
is the SO(n) AKLT model 0, [13, Isll, whose ground 
state can be represented as a projected BCS state, by 



replacing gij = 1/ (zi — Zj) in ([T]) with gij — 1 1^. For 
the Ising disordered phase, the ground state of the spin 
chain is dimerized [20j and hence the valence bonds are 
short-ranged. In this case, a proper Cooper pair wave 
function for the projected BCS state can be chosen as 
gij ~ exp(— |zi — Zj\/£,), where ^ is the length scale of the 
valence bonds. In the extreme Cooper pair wave 

function that is nonvanishing only between neighbor- 
ing sites yields Majumdar-Ghosh-like state, which cor- 
responds to perfect dimerization. These results imply 
that both Ising ordered and disordered phases close to 
SO(n)i criticality are well described by projected BCS 
states with properly chosen g^j. Indeed, for n = 3, it was 
shown [13] that the projected BCS states with Cooper 
pair wave functions generated from Kitaev's Majorana 
chains [s^l are good variational wave functions for the 
Haldane (Ising ordered) and the dimerized (Ising disor- 
dered) phases. 

2D chiral spin liquids.- After establishing the relevance 
of projected BCS states H]) for SO(n)i criticality in ID, 
we move on and discuss their properties in a 2D square 
lattice, where the z's in ^ are now complex coordinates 
of the lattice sites. In an analogy with fractional quan- 
tum Hall (FQH) states constructed by conformal blocks 
of their gapless edge CFTs 33, [s^, the chiral correla- 
tor Q from SO(n)i WZW model (n massless Majorana 
fermions) yields chiral spin liquids, which break time re- 
versal symmetry and are spin counterparts of FQH states 
[s^]. From the projected BCS form ([1]), the Cooper pair 
wave function l/{zi — Zj) now corresponds to the topolog- 
ical phase ofp+ip superconductors [36| supporting chiral 
gapless Majorana edge modes, which justifies the above 
bulk-edge correspondence. Below we focus on the any- 



onic quasiparticle excitations in these 2D states, which 
have intriguing properties depending on the parity of n. 
A more detailed analysis of these 2D states will be pre- 
sented elsewhere. 

Interestingly, the quasiparticles built upon the SO(n) 
states support non-Abelian statistics for odd n. Let us 
adapt the CFT approach of creating quasihole excitations 
in FQH states [33| to our spin system. For odd n, the 
SO(n)i WZW model has three primary fields: identity 
field / (singlet), vector field v (Majorana fermions), and 
spinor field s. Following the CFT approach, creating 
quasiparticles in the SO(n) state is achieved by adding 
spinor fields s in the chiral correlator Then, the 

statistics of quasiparticles are encoded in the fusion rules 
of the primary fields. In fact, the spinor fields have a 
nontrivial fusion rule sxs = I+v, together with sxv = s 
and V X V = I. These fusion rules resemble those in Ising 
CFT ((TX(T = / + e, (Txe = cr, and e x e ~ I), which are 
resp onsible for the non-Abelian statistics of Ising anyons 
[37| . Indeed, the Majorana free field representation of 
SO(n)i WZW model allows us to identify the spinor fields 
s with conformal weight hg = n/16 as a product of n 
Ising a fields {he — 1/16). Thus, we conclude that the 
SO(n) states support non-Abelian Ising anyons for odd 
n. Note that the case with n = 3 recovers the physics of 
the Moore- Read states H, 11] , while for odd n > 5 they 
are natural generalizations of the Moore- Read states. 

Now we show that the SO(n) states only support 
Abelian anyons for even n. This subtle difference roots 
in the fusion rules of the SO(n)i primary fields. In con- 
trast to odd n case, the SO(n)i WZW model with even n 
has two spinor primary fields s+ and s_ [isl ]. apart from 
the usual identity and vector fields. The fusion rules of 



spinor and vector fields are 



S- and S- X V = s+. 



Depending on the parity of n/2, the fusion rules involving 
two spinor fields are s+ x s+ = s_ x s_ = /, s+ x s_ = w 
for even n/2 and s+ x s+ — x — v, s+ x — I 
for odd n/2. However, due to the absence of multiplicity 
in the fusion outcome, only Abelian anyons can exist in 
the SO(n) states with even n. 

Conclusion and perspective.- To conclude, we have 
proposed a class of projected BCS states and derived 
their parent spin Hamiltonians. These wave functions 
also have an infinite MPS form generated by chiral cor- 
relators of Majorana fermions. In ID, they can be viewed 
as SO(n) generalizations of Haldane-Shastry models and 
capture the physics of SO(n)i WZW model. These re- 
sults indicate that modified projected BCS states are 
good variational ansatz for describing Ising ordered and 
disordered phases close to SO(n)i criticality. In 2D, the 
SO(n) states are chiral spin liquid states, which support 
non-Abelian Ising anyons for odd n and Abelian anyons 
for even n. An open question that deserves further inves- 
tigation is whether these 2D chiral spin liquids are rel- 
evant for physical models and materials [39[. Moreover, 
our 2D toy models may also shed light on an important 
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open issue in condensed matter: Can superconductivity 
arise after doping these antiferromagnets? 
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Supplemental material 

Equivalence of the projected BCS state and the infinite MPS 

In this Section, we prove the equivalence of the projected BCS state and the infinite MPS. 

Projected BCS state 

Let us first expand the projected BCS state 

|0) 



\i<j * ^ a=l y 



= PcU 



N 



a=l 



^-^ ^-^ Zi — Zi x\\a 

Na=0 {Na even)^(<.)<...<^(a) 



,t 

Zj' x{~',a ^n''«- 



'•C (a) 



|0) 



E E E ••• E npf«(7:^) 

X(c^,> •••c'^rii )(c^2i ■■■c'(2) )---(c^„) •••c''',„N )|0) 



where Na {a = 1, ... ,n) is the number of fermions (i.e. |n") vector in the configuration), a;^"^ < • • • < are 
the positions of fermions in the lattice, and the Pfafiian factor PfaC ^.l^^. ) is restricted to the positions for the 
fermions. 

The next step is to implement the Gutzwiller projection. Note that the Gutzwiller projector Pq requires single 

occupancy. Thus, the positions of fermions (x^j^'' < ■ ■ ■ < x^j, x^2 ^ < ■ • • < x^^, ... , .tJ"-* < • • • < a;^'*) must be 
all different from each other, so that each site has exactly one fermion. As a result, we have J22=i -^a ~ After 
implementing the Gutzwiller projector, we obtain 



i*)= E E npM^) 

JVi,iV2,...,iV„=0 iNi+N2 + -+N^=N and oven) allowed a:*"' <---<2;<^_^ "=1 * ^ 

X(c\i) )(C^(2) ■••C^a) )---{c\n) ■■■c'm )|0) 

The final step is to rearrange the positions of fcrmionic operators so that they can be identified as a spin state. 
This rearrangement only results in a sign, depending on the positions of fermions 



w= E E npu^) 

JVi,JV2,...,JV„=0 (JVi+JV2 + -+JV„=JV and Na even) allowed x'f^ <-<x^Sl ' ^ 

xse-nrr^^^ ^(1) ^(2) ^(2) r^^'^Mr^^) r'^^ r^^^ r'^^ r^") r^"^^ 

'^°B"V'''l ) • • • 5 -''jVi ' 1 5 ■ • • 5 •''iV2 ' ■ ■ ■ ' 1 ' ■ ■ ■ ' •''Ar„ / 1 1 ' ■ • • I •*'A''i ' '''1 ) • • ■ I •''JV2 ' ■ ■ ■ ' 1 ' ■ ■ ■ ' JV„/ 

where jxj^' , . . . , x^j , x^^^ , . . . , x^^ , . . . , xj^' , . . . , x^^) is a spin configuration labeled by the positions of the vector jn") 

and sgn(x^^\ . . . , xj^\ . . . , xj^\ . . . , x^"\ . . . , xj^^) is the signature of permutation due to the sign factor coming 

from fermionic anticommutation relations. 

Thus, the projected BCS wave function can be written as 

n ^ 

*({xW},{x(2)},...{x(")}) = sgn(xl^\...,x(^^xf\...,x(^],...,xl"\...,x(^))n^^^ (1) 

a=l 
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where {x^""^} is the set of positions satisfying x'f'' < • • • < xj^'' {Na even and X]a=i = N). 

Infinite MPS 

Let us now consider the infinite MPS 

n 

ai .. . . ,a7V" — 1 

where ^(ai, . . . , un) ai'e given by the chiral correlators of Majorana fermion fields (a = 1, . . . , n) 

vfCai, ...,aN) = (x"^ (^i)x"^ (^2) • • • X"" (^Jv)) 
To evaluate ^(ai, . . . , oat), we use the two-point correlator of Majorana fermions 

z — w 

The multipoint correlators of Majorana fermions are obtained by Wick's theorem 



(X'^(21)X'^(^2)---X"(^A^J> = 



Pf«(3-^) ^aeven 
Na odd 



Therefore, to obtain a nonvanishing ^'(ai, . . . ,aAr), we must have even N. Additionally, Na, the number of vectors 
|n") in the spin configuration, must also be even for all a = 1, . . . , n. 

To compare with the projected BCS wave function, let us evaluate the superposition coefficient of the infinite MPS 
for a spin configuration, which has Na vector |n") at positions a;^°^ < ••• < a;^^ {Na even and J22=i = N). 
Taking into account the anticommuting nature of Majorana fermion fields, we first pick up the vectors |n^) in the 
spin configuration and rewrite the infinite MPS as 

{^(2)}, . . . = sgn(xW, . . . , a^j^j , yi, . . . , yA^-L J(x"=^(^,(i,)x"='(^4i)) ' • ' x"='(^,a) )) 

x(/(%Jx''(%.) (Ml) 

/ (1) (2) \-n^ r ^ 

= sgn{x\ ' , . . . ,xy^,yi, . . . ,yN-Li)Pia=i[ — - — ) 

Zi Zj 

where the positions yi < ■ ■ ■ < yN-Ni correspond to the vectors |n^) with b ^ 1. The above steps can be repeated 
from b = 2 to n. In the end, we obtain 

n 

*({xW},{x(2)},...{a:(")})=sgn(x(^\...,xW,4^\...,xi^],...,x("\...,4")^ (2) 

T Zi Zj 

Comparing with Eq. (H]), we conclude that the infinite MPS and the projected BCS state are equivalent. 

Derivation of tlie parent Hamiltonian 

In this Section, we derive the parent Hamiltonian for the infinite MPS. 

Brief summary of tlie SO(n)i WZW model 

For infinite MPS associated to WZW models, the derivation of the parent Hamiltonian relies on the existence of 
null vectors in the representation spaces of Kac- Moody algebra [9]. For SO(ri)i WZW model, the Kac-Moody algebra 
is defined by 

[Jf, ] = ^r^'^^'^-^^^i™ + ndab,cdSn+m,o n,meZ (3) 
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where repeated indices are summed over and the SO(n) structure constant f^b.cd.ef gjygj^ i-,y 

jab,cd,ef ^ SadSbeScf + hc^ae^df " SacSbe^df - SbdSae^cf 

For odd n (n > 3), the SO (71)1 WZW model has three primary fields respectively in singlet (denoted by /), vector 
(v) and spinor representation (s), whose conformal weights are hi = 0, hy = 1/2 and hs = n/16, respectively. For 
even n (n > 4), apart from the primary fields in singlet and vector representations {hj = and hy = 1/2), the SO(n)i 
WZW model has two primary fields in spinor representations (denoted by s+ and s_), whose conformal weights are 
hs+ = hs_ — n/16. The SO(n)i WZW model has central charge c — n/2 and can be constructed by combining n 
Ising models (c — n x i). 

For both odd and even n, the primary fields in the vector representation have conformal weight hy = 1/2 and are 
naturally interpreted as Majorana fermions, which are the key ingredients for us to construct the infinite MPS. 
For each Majorana fermion (a = 1, . . . , n), a primary state |x") can be defined by 

lx'^) = x'^(o)|o) 

where |0) is the vacuum of the WZW model and satisfies J,°>o|0) = 0. When acting on the Kac- Moody currents, the 
primary states satisfy 

n 
d=l 

Jn'lxl = (n > 0) (4) 

where L"'' are given by 

[L"'^)cd = i{SacSbd - SbcSad) 

Note that L"'' form a closed SO(n) algebra 

_ ^rab,cd,ef j^ef 



Null vectors and parent Hamiltonian 

To derive the parent Hamiltonian, one has to find the null vectors in the SO(n)i Kac-Moody algebra. The null 
vectors are descendant states satisfying 

J,f|</.)=0 (n>0) 
For our purpose, we look for null vectors with the following form; 

a<.h,c 

where W^^^ are the coefRcients that have to be determined. They satisfy the orthonormal condition 

E (WabcTWtbc = Sa'd 
a<b^c 

In general, the tensor Wq^c corresponds to a Clebsch-Gordan decomposition. For the SU(2)fe WZW model, the SU(2) 
Clebsch-Gordan coefficients are known Q. However, we are not aware of a closed form for the SO(n) Clebsch-Gordan 
coefficients. To overcome this difficulty, let us consider the norm of 

i^'i^") = E T.(^a'b'c'rwUx''\jf''r\\x'') 

a' <b' ,c' a<b,c 

= E E WUc'Ma'b'CabcW^bc 
a' <b\c' a<b,c 
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where is viewed as a column vector and M is a matrix defined by 

Ma'b'c'.abc={x''\jf''r\\xl 

If 10'') is a null state, | = {W^y MW^ — 0. Since M comes from the norm of two descendent states, it is 
a positive-semidefinite matrix satisfying {W^y AIW^ > 0. Therefore, identifying the orthonormal vectors that 
belong to the kernel of M gives us all null states 10''). For our SO(n)i WZW model, we can write down the explicit 
form of M 

= ix^'iiJi'^^j-'im 



e/ g 
ef 

where we used Kac- Moody algebra and the properties of the primary state (|H). In this way, the null vectors for 
the SO(n)i WZW model are obtained. 

Let us mention that the above approach is a systematic way of finding null vectors and can be easily generalized 
to other WZW models. The role of the positive-semidefinite matrix M is similar to the Gram matrix for defining the 
Kac determinant [l5l | in conformal field theory. 

After obtaining W^^^ for all the null vectors, we define the following K tensor Q 

K'bV' - E«'&'c')*W^at (« < b and a' < 5') 

d 

Let us write these K tensors as n x n matrices, K'^^^'i ^, ~ (i^°/'b/)c,c'- Then, K'^';'^^, have the following compact form 



T^ab ^ c Tl ^ f ab ,a' b' ^cd T cd i ^ ^ ( y ab y a'b' i j a' b' j ab\ 

3 6(n — 1) 6(n — 1) 



Following Ref. we define an operator Kf' (\ < a <h <n) 



N 



] = l{^i) c<d 
N 



E 



N 
E 



-if - 



n + 2 



3 ^ 6(n-l) 



■j^ab,cd,ef^f ^cd ^ 



( T ab T cd I T cd T ab\ J- cd 



3 ^ n-l 



Lf{L,-L,) + -{L,-L,)Lf 



where = {zi + Zj)/{zi ~ Zj) and L,; • = J2a<b^i^^'j^- These operators annihilate the infinite MPS, Af jvf) = 
Vi, a, 6. Moreover, the infinite MPS is an SO(n) singlet and therefore J^i 1^) = ^o, b. Then, an SO(ri) symmetric 
parent Hamiltonian can be defined by 

H=Y1 (Af )^Af + ^E(E ^f)' +^0 ( J > 0) 

i,a<b a<b i 

whose ground state is the infinite MPS with energy £"0. 
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In ID, we use Zj = exp(i^j) to ensure translational invariance. Choosing J = 2(A'' — 2)/3 and Eq = 
l)iV(iV^ — 4)/9, we arrive at the following explicit form of H: 

H = -y^ w^j [^^^ {Li ■ Lj) + r,? \ \ • -^j)^] ~ qII \\ WijWik{Li- Lj){Li- Lk). 

To obtain the above form, the following identities are quite useful 

Y,Lf{U.L,)Lf = U-L, {i^j) 

a<b 

J2 Lf>{U ■ hfLf = 2(n - 1) - (n - 2)(L, • Lj) - {U ■ Ljf {i + j) 

a<b 

^ Lf{Li ■ Lj){Li ■ Lk)Lf = 2{Lj ■ 4) - {U ■ Lk){Li ■ L,) (i ^ j ^ k) 



-2{n - 
(5) 



a<b 



Jastrow and PfctfHan wave functions in Cartan basis 

In this Section, we derive the explicit Jastrow and Pfafhan forms of the wave functions in Cartan basis. 



Cartan basis 

Let us first define the Cartan basis. 
The SO(n) algebra is defined by 

For n = 21 and 2/ + 1, we can choose at most / (rank of the algebra) mutually commuting generators as 
...,L2;-i.2i. In the vector basis, the SO(n) generators are defined by L"^ = i(|n")(n''| - \n''){n''\) 
{1 < a < b < n). Diagonalizing the Cartan generators gives us the following Cartan basis: 



|1,0,... 
1-1,0,... 
|0,1,0,... 
|0,-1,0,... 

|0,0,...,1 
|0,0,...,-1 



V2 

1 

71 

1 

1 

71 

1 

71 
1 

71 



For S0(2/), the above basis is already complete. For S0(2/ + 1), we have an additional vector |n^'+^), which is 
annihilated by all Cartan generators. Thus, we have the following extra vector for S0(2Z + 1): 

|0,0,...,0) = |n2'+i) 

Thus, the Cartan basis for S0(2Z) and S0(2Z + 1) can be compactly written as 

|0,...,m„ = ±l,...,0) = -^(|n2«)±i|n2"-i)) {a = l,...,l) 

|0,0,...,0) = In^'+i) 
Note that iria is the eigenvalue of the Cartan generator i2a-i,2a_ 
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Even n = 2l 



Now we derive the Jastrow and Pfaffian forms of the SO(n) wave functions in the Cartan basis. Actuahy, this goal 
can be achieved from either the projected BCS form or the infinite MPS form. In the following we use the projected 
BCS form to derive the results. 

For S0(2I), we define the slave-fermion operators in the Cartan basis 



^mi=±l ~ 



Cm2=±l 



V2 



(4( ^ *4i-i) 



After changing the basis, the S0(2Z) valence-bond singlet operator in the projected BCS state is written as 

21 I 

EJ J _ sr^fJ J ,J J \ 

a=l a=l 

Using the above form, the S0(2Z) projected BCS state is rewritten as 
l*so(20) = ^Gexp(5:-^f:4ctJ|0) 

\i<j ' ^ a=l J 

= Pg exp 4,m^ = l4,m^ = -l |0> 

a=l / 



n 

I 



~c c 

. i^ntoc — l j,mot 



1 |o) 



= PGn 



N/2 

E E 



V detf^— ) 



X c 



is the determinant of the x Cauchy matrix restricted to the positions of 



• • • c 



„(») 



|0) 



where det 



(pr-pW"i).(9i"'-4"i) 



4i„=i ^iid '^mo,=-i fermions. The following useful identity reduces the Cauchy determinant to a product of Jastrow 
factors: 



det 



1 



(_l)^JV„(iV„-l) 



Note that the sign factor in the Cauchy determinant can be absorbed by rearranging the fermionic operators 



(_l)iiV<,(iV„-l)(^t ,t ...,t )(,t ,t ...,t ) 

p ' mo = l pi ^m<:, = l B*°\mc = l or',Tnc = -l oi ',mc = -l oi°' mo = -l 



„(a) 



9i 



92 
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Therefore, we obtain 



I* 



SO(2/)/ 



Pen 



N/2 

E E 



E 



Pi 



X (c^ 



|0) 



In the next step, we collect the positions p " < ■ ■ ■ < and ql"' < ■ ■ ■ < gj^j into a single set {a;^")} with 
a;^"'* < • • • < Xj"^ . Then, the Jastrow factors can be written as 



- n 

l<i<i<2Arc, 



up to a sign factor. However, the sign factor can again be compensated by rearranging the fermionic operators to the 
correct order according to x'f ^ < ■ ■ ■ < x'^ . 

The last step is to implement the Gutzwiller projection and switch to the spin basis. As a result, we obtain the 
S0(2^) wave function in Cartan basis 



*(w)=prnrin(^ 



(6) 



a—1 i<j 



where pm — sgn(a;$^"'^'' , . . . , , • • • , x^' , . . . , x'"^' ) if J^i ^a.i = Va and pm = otherwise. 

For n = 4, one can further simplify Eq. and show that it is equivalent to a product of two spin-1/2 Haldane- 
Shastry states, if the four vectors are interpreted as two spin-1/2 states. 

Odd n = 2l + l 

Comparing to S0(2/), we have an additional slave-fermion operators in the Cartan basis for S0(2Z + 1) 

Then, the valence-bond operator for S0(2I -I- 1) is expressed as 
21+1 I 

0=1 Q — 1 

The expansion of the S0(2Z + 1) projected BCS state is given by 

/ 21 + 1 \ 

l*SO(2^+i)) = Pcexp -^-T^. E 4aS> |0) 

\i<] ' J a=l j 

= Pg exp -^-^cJ_„^oc] E T~r E Am^=A.m^=^l |0) 



N 



No=0 (No even) ^(0)<... 2,0') 



I (0) n (0) n 



( E E 4™„=icL„=-i ) 

* a=l j 



|0) 
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where the positions of the extra fermion c„j^q are labeled hy x\ ' < ■ ■ ■ < x)^^. The rest of the calculation is similar 
to the S0(2/) case, except for the presence of a Pfaffian factor due to the extra fermionic mode cJ^^q- After some 
algebra, we obtain the S0(2Z + 1) wave function in Cartan basis 

1 ' 

*({m}) = p,„Pfo(^— ) H - z,)"-*""- (7) 

a— 1 i<j 

where pm = sgn(a;^°^ . . . x^^^ , x[^^ , . . . , z^j , . . . , x^'^ ,x^^^) if J2i = Va and pm = otherwise. 

For n = 3, the signature of the permutation in pm is reduced to the Marshall sign in the spin-1 Haldane-Shastry 
state [9|]. 



